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1 Introduction 



The standard field theoretical perturbative approach to quantum gravity (see 
the introduction to this subject in [jlj and references therein), considered as a flat 
space-time relativistic quantum field theory of gravitons, massless rank-2 tensor 
fields, coupled to photons, massless vector fields, leads to non-renormalizable 
ultraviolet (UV) divergences. These were found by means of dimensional regu- 
larization and background field method 0, ||, ||] and ||. 

A closer investigation of loop calculations in quantum field theory shows that 
the reason for the UV divergences lies basically in the fact that one performs 
mathematically ill-defined operations, when using Feynman rules for closed loop 
graphs, because one multiplies Feynman propagators as if they were ordinary 
functions. 

We show how it is possible to overcome these discouraging outcomes by 
applying an improved perturbation scheme which has as central objects the 
time-ordered products and as constructing principle causality. The S-matrix is 
constructed inductively as a sum of smeared operator-valued n-point distribu- 
tions avoiding UV divergences. 

UV finiteness of S'-matrix elements is then a consequence of a deeper math- 
ematical understanding of how loop graph contributions have to be calculated. 

This idea goes back to Stuckelberg, Bogoliubov and Shirkov and the pro- 
gram (causal perturbation theory) was carried out successfully by Epstein and 
Glaser || for scalar field theories and subsequently applied to QED by Scharf Q , 
to non-Abelian gauge theories by Diitsch et al. || and to quantum gravity ||. 

In this work we carry out the analysis of the coupled quantized Einstein- 
Maxwell system, Sec. 2.2. Two main topics will be investigated in this paper: 
the gauge structure of the theory and the UV finiteness of loop graphs in second 
order perturbation theory. 

The first is formulated by means of a 'gauge charge' that generates infinites- 
imal gauge variations of the fundamental free quantum fields. Gauge invariance 
of the S'-matrix implies then a set of identities between the C-number part of 
the 2-point distributions, Sec. 2.4, which implies the gravitational Slavnov-Ward 
identities (SWI) 0. 

The second is obtained by applying the 'causal perturbation' scheme, Sec. H 
to the calculations of the photon loop graviton self-energy, Sec. and of the 
photon self-energy, Sec. ||, in second order perturbation theory. In both cases the 
results are UV finite and cutoff-free and the inherent ambiguity in the normaliza- 
tion of the 2-point distribution is settled by requiring appropriate normalization 
conditions for mass and coupling constant. 

In addition, the causal scheme preserves the gauge symmetries of the theory. 
The trace of the graviton self-energy tensor vanishes. This property was mani- 
festly broken using dimensional regularization in order to extract finite results 
from UV divergent quantities |ll]] . 

From the photon loop graviton self-energy we sketch the derivation of quan- 
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turn corrections to the Newtonian potential between two spinless massive bodies 
in the static non-relativistic limit, Sec. fl5| . 

Gauge structure and one-loop calculations with gravitational self-coupling 



are not considered here, see |9|, |l2j, [|13| , and [|14|. These references contain 



also the notations and the conventions used here. 

We use the unit convention: h = c = 1, Greek indices a, /3, . . . run from to 
3, whereas Latin indices i,j, . . . run from 1 to 3. 

2 Quantized Einstein— Maxwell System and Pertur- 
bative Gauge Invariance 

2.1 Inductive S-Matrix Construction 



In causal perturbation theory py, [15|, one makes an ansatz for the 5-matrix 
as a formal power series in the coupling constant, namely as a sum of smeared 
operator-valued distributions: 

°° 1 f 

S(g) = 1 + ^ — \d A xi...d A x n T n (xi,... ,x n )g(x{) ■ g(x n ) , (2.1) 

n=l 



where the Schwartz test function g 6 5(M 4 ) switches adiabatically the interac- 
tion and provides a natural infrared cutoff in the long range part of the interac- 
tion. The S-matrix maps the asymptotically incoming free fields on the outgoing 
ones and it is possible to express the T n 's by means of free fields. Interacting 
quantum fields are never used in the causal scheme. 

The n-point distribution T n is a well-defined 'renormalized' time-ordered 
product expressed in terms of Wick monomials of free fields : 0{x\, ... , x n ): 

T n (%\, . . . , Xn) — ^ '. 0{xi , . . . , Xn) . t n ( y X\ X n , • • • , X n —\ Xn) . (2.2) 
O 

The tnS are C-number distributions. T n is constructed inductively from the 
first order Ti(x), which plays the role of the interaction Lagrangian in terms of 
free fields, and from the lower orders Tj, j = 2, . . . , n — 1 by means of Poincare 
covariance and causality. The latter leads directly to UV finite and cutoff-free 
redistributions . 

2.2 Quantized Einstein— Maxwell System 



In the context of linearized gravity (without graviton self- interactions ||, [ 13 1 



[14 1) coupled to photon field Q, the interaction between the quantized symmet- 
ric tensor field h a ^(x), the graviton, and the quantized vector field A^{x), the 
photon, is considered in the background of Minkowski space-time. 
The free graviton field satisfies the wave equation 

ah^(x) = 0, (2.3) 
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and is quantized according to 

[h^(x),h a ^(y)] = -ib^D (x-y), (2.4) 
where the fr-tensor is constructed from the Minkowski metric 

y*>*f> ._ 1 (jfajfP + ^rf" _ rj^rj ^ , (2.5) 

and Dq(x) is the mass-zero causal Jordan-Pauli distribution: 



D (x) = D ( +\x) + D l Q -\x) = ^S(x 2 ) sgn(x°) 



(2vr) ; 



■ jd 4 p5{p 2 )sgn(p )e-~ ip:1 



(2.6) 



(see [12] and [13] for the details of this procedure). The additional degrees 
of freedom present in the symmetric tensor field (gravitons are massless 
spin-2 particles) could be eliminated by imposing appropriate gauge and trace 
conditions, but these are disregarded in the following and only considered later 
as a characterization of the physical states |l^] . 

From the point of view of Lagrangian field theory |J, the Hilbert-Einstein 
Lagrangian C HE = —2 k~ 2 yj—g g^ u R^u, (k 2 = 32irG), is written by means of 
the Goldberg variable g^ v := yj— g g^ v and expanded as g^ v = rf- v + kW 
around the flat space-time rf v = diag(+l, —1, —1, —1). From the zeroth order, 
choosing the gauge = 0, one obtains (^3]). 

The free photon field fulfils also the wave equation 

nA»(x) = Q, (2.7) 

and is quantized as 

[A^{x),A v {y)] = +i V ^D (x-y). (2.8) 
Also the field strength F^ u := d^A v — d v A^ will be used in the following. 



2.3 Gauge Invariance and First Order Coupling T^(x) 

The 17(1) em gauge content j/J of the photon field is formulated by means of the 
gauge charge Q A 

Q A := f d 3 xA v (x), u d° v(x), (2.9) 

x°=const 

where v{x) is a C-number scalar field satisfying Ov(x) = 0. Actually, due to the 
generality of the gravitational coupling, one should also consider the coupling 
between h^ v and the quantized electromagnetic ghost v and anti-ghost v scalar 
fields, but this would have the form Kr) fiv h^ 1 ' vdv from \/—g vOv and therefore 
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can be written as a total divergence in the sense of vector analysis because of 
the presence of two equal derivatives on three massless fields. For this reason 
the electromagnetic ghosts will not be considered here. 

On the other side, the gauge content of the graviton field (which is related 
to the general covariance of g^ u (x) under coordinate transformations |16|] ) is 
formulated by the gauge charge Q, |l4j 

Q:= J d 3 xh >lu (x) !U d^ u M (x). (2.10) 

x°=const 

In order for gauge invariance to first order in pure quantum gravity to hold 



(see Eq. ( 2.15| )) |§, for the construction of the physical subspace [12] and in 



order to prove unitarity of the 5-matrix on the physical subspace fl2|| , we need 
to quantize the ghost field u^^x), together with the anti-ghost field u u (x), as 
fermionic vector fields 

□u^(x) = 0, Du v (x)=0, {u»{x),u v {y)} = iff D Q {x - y) , (2.11) 

whereas all other anti-commutators vanish. The gauge charges generate the 
following infinitesimal gauge variations of the fundamental quantum fields: 

dQhT{x) := [Q,h^(x)] = -ib^uP{x) ta , 
d Q u a (x) := {Q,u a (x)} =0, d Q u a (x) := {Q, u a (x)} = i h a/3 (x) i/3 , 
d QA A a (x) := [Q A ,A a (x)] =id%v(x), d QA F a ?(x) = 0. (2.12) 
The operator obeys also the Leibniz rule 

d Q (AB) = (d Q A)B + (-l) n ^Ad Q B, (2.13) 

for arbitrary operators A and B, where n G {A) is the number of ghost fields 
minus the number of anti-ghost fields in the Wick monomial A. 

For convenience of notation, the trace of the graviton field is written as h = 
h?~/ and all Lorentz indices of the fields are written as superscripts whereas the 
derivatives acting on the fields are written as subscripts. All indices occurring 
twice are contracted by the Minkowski metric rj^ . We skip the space-time 
dependence if the meaning is clear. 

The first order interaction between graviton field and photon fields is dictated 
by gauge invariance and by the assumption that the graviton interacts with the 
traceless and conserved electromagnetic energy-momentum tensor, see below. 

Gauge invariance of the 5-matrix means formally 

lim dnS(g) = 0, and lim dn A S(g) = . (2.14) 

Since the existence of the adiabatic limit g — > 1 in massless theories can be 
problematic, we consider the condition of perturbative gauge invariance to n-th 
order 

dQT n (xi, ... , x n ) = dQ A T n (xi,. . . , x n ) = sum of divergences , (2.15) 
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which implies S-matrix gauge invariance, because divergences do not contribute 
in the adiabatic limit g — > 1 due to partial integration and Gauss' theorem. 
We define T^(x) as 

Tf(x) :=t£ :h^{x)b, vap Tf M {x):, (2.16) 

where T FM { X ) corresponds to the electromagnetic energy- momentum tensor and 
reads 

Tf M {x) := -F 0f r(x)F^(x) + \<f# F>»(x) F pa (x) . (2.17) 

Actually, the 6-tensor is not essential, because T FM is traceless. Its origin lies in 
the Goldberg variable expansion. T±(x) corresponds to the term of order k in 
the expansion of the generally covariant electromagnetic Lagrangian density 

C A = -\^F^F aP g^g^ , (2.18) 

as a function of g^ u around the flat space-time. 

Perturbative gauge invariance to first order is readily established. T*{x) is 
U{\) EM gauge invariant 

d QA T 1 A (x) = 0, (2.19) 

due to d QA F^(x) = 0. The 17(1) em gauge structure of the theory is rather 
trivial and will not be considered. 

Gauge invariance with respect to the gauge charge Q 

d Q T?{x) = I :u"{x), a T^ M {x):= :u'(x)T&:) =: d^T^x) (2.20) 

holds because of T"£(x) 5j g = 0. T°, x (x) is the so-called Q-vertex Q, ||. It 
allows us to formulate the condition of perturbative gauge invariance to the 
n-th order in a precise way : 

n d 

dQT n (x\,... ,x n ) = 22 q~i; T n/i( x i> ■ ■ ■ > x z>--- > x n), (2-21) 
i=i X[ 

where T%n is the 'renormalized' time-ordered product, obtained according to the 
inductive causal scheme, with one Q-vertex at xi, while all other n — 1 vertices 
are ordinary T\ -vertices. 

2.4 Consequences of Perturbative Gauge Invariance to Second 
Order 

Without performing any calculation, from the structure of we can anticipate 
that the two-point distribution describing loop graphs has the form 

T 2 (x,y) l0 °P s =:h^(x)h^(y): t hh {x,y) a ^ v + :F^{x)F^{y) : t FF {x,y) a ^ v . 
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Here, the first term represents the photon loop graviton self-energy and the 
second term the photon self-energy. The C-number distributions thh an d t FF 
will be explicitly calculated in Sec. |I| and in Sec. [|, respectively. 

Perturbative gauge invariance to second order, Eq. ( 2.21 ) with n = 2, enables 
us to derive a set of identities for these distributions by comparing the distri- 
butions attached to the same external operators on both sides of Eq. ( |2.21| ), as 
in 0. Calculating d Q T 2 (x , y) loops we obtain 



d Q T 2 (x,y) loops =+ :u p (x) :<7 h^(y): ( - i b a ^ t hh (x, y) aPlu/ ) + 
+ : h^(x)u p (y), a b pvp ° t hh (x, y) a ^ v ) . 



(2.22) 



On the other side, the loop contributions coming from T£^(x,y) can only be of 
the form 



T° n (x,y) =:u p (x)h^(y): t^x, y) m „+ :u° (x)h p » (y) : t uh (x,y)^ . 

Applying d% to the expression above we find 

ggTZ/^y) =+ :u p {x)^h pv {y): {tZ h (x, y) PVbV + r, p a t uh (x, y)^} + 
+ :u p {x)h pv {y): 8 x a {t a uh {x, y) pia/ + rj P t uh {x } y)„ v } . 



(2.23) 



(2.24) 



We compare the C-number distributions in ( |2.22| ) and in ( 2.24 ) attached to the 
external operators 

:u p {x) >a h^{y): and :u p {x)h pv (y): . (2.25) 
No such terms come from T 2 , 2 (x,y). Therefore, we obtain the identities 

-i b p ™P t hh (x, y) a ^ v = {t° uh (x, y)% + rT y)^} , 

= % {t° lh (x,yy ixu + ri< } °t uh (x,y) tlu } . 

By applying d x to the first identity and inserting the second one, we obtain 

b a ^ pa d x a t hh (x,y) a ^ u = (2.27) 

This identity has been explicitly checked in Eq. ( |4.20D and implies the grav- 
itational Slavnov-Ward identities Q, [jlO| for the two-point connected Green 
function (see Sec. 4.3) 



3 Causal Construction of the Two-Point Distribution 

We outline briefly the main steps in the inductive causal construction of T 2 (x, y) 
from the first order interaction T A (x). Following the inductive scheme we 
first calculate 

R' 2 (x,y) := -T 1 A (y)T 1 A (x), A' 2 (x,y) := -Tf(x) Tf(y) . (3.1) 
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From these two quantities, we form the causal distribution 



D 2 (x,y) := R' 2 (x,y) - A' 2 (x,y) = [T 1 A (x),T 1 A (y)} . (3.2) 

In order to obtain D 2 (x, y), one has to carry out all possible contractions between 
the normally ordered T\ using Wick's lemma, so that D 2 {x, y) has the following 
structure (see 

D 2 (x,y) = Y J -0(x,y): d%(x-y). (3.3) 

o 

d 2 (x — y) is a numerical distributions that depends only on the relative coordi- 
nate x — y, because of translation invariance. 

D 2 {x,y) contains tree (one contraction), loop (two contractions) and vac- 
uum graph (three contractions) contributions. Due to the presence of normal 
ordering, tadpole diagrams do not appear. D 2 (x, y) is causal, i.e. supp(c?2 (z)) Q 
V+(z) U V~(z), with z := x-y. 

In order to obtain T 2 (x,y) we have to split D 2 (x,y) into a retarded part, 
R2(x,y), and an advanced part, A 2 (x,y), with respect to the coincident point 
z = 0, so that sup-p(R 2 (z)) C V+(z) and supp(A 2 (z)) C V~{z). This split- 
ting, or 'time-ordering', has to be carried out in the distributional sense with 
great care so that the retarded and advanced part are well-defined and UV 
finite §, §. 

The splitting affects only the numerical distribution d 2 (x — y) and must be 
accomplished according to the correct 'singular order' uj(d 2 ) which describes 
roughly speaking the behaviour of d 2 (x — y) near x — y = 0, or that of d 2 (p) in 
the limit p — ► oo, respectively. If to < 0, then the splitting is trivial and agrees 
with the standard time-ordering. If a; > 0, then the splitting is non-trivial and 
non-unique: 

<"(d?) 

d°(x-y)^r°{x-y)+ £ C a>0 D a 5™ (x - y) , (3.4) 

|a|=0 

and a retarded part r 2 {x — y) is obtained in momentum space by means of a 
dispersion-like integral of the type 

which requires modifications in the case of massless theories B. 



Eq. (3.4) contains a local normalization ambiguity: the C aj0 's are undeter- 
mined finite normalization constants, which multiply terms with point support 
D a <5 (4) (x _ ^jja j g a p ar ^ a i differential operator). This freedom in the nor- 
malization has to be restricted by physical conditions. For example, Lorentz 
covariance and gauge invariance will be used in our discussions. 
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Finally, T 2 {x,y) is obtained by subtracting R' 2 (x,y) from R 2 (x,y) and can 
be written as 



T 2 (x,y)+iV 2 (x,y) = ]r:0(x,y): i #(s - y) + ^ ^ £> a 5 (4) (x - y) , 

^ |a|=0 ' 



(3.6) 



where N2(x,y) represents the normalization terms given by the second term on 
the right side. 

4 Photon One-Loop Contribution to the Graviton 
Self-Energy 

4.1 Causal D 2 (x, y) -Distribution for the Photon Loop 

We now apply the causal scheme described in Sec. |3| in order to calculate the 
2-point distribution T 2 (x,y) sSE that describes the photon loop graviton self- 
energy. 



First of all, from the commutation rules (2.4) and ( |2.8| ) we derive the con- 
tractions between two field operators 

C{h aP {x) h^(y)} := [h aP {x) { '\h^ u [y) m ] = -ib a ^ u D { +\x - y) , 
C{A^x) A v {y)} := [A» (x)", A" (y)™] = +i rT D«\x - y) , (4.1) 

where (±) refers to the positive/negative frequency part of the corresponding 
quantity. 

The distributions R' 2 (x,y) 9SE and A' 2 (x,y) 9SE , defined in Eq. (|3~T|) , are ob- 
tained by performing two photon contractions between the two first order inter- 
actions Ti(x) and T^{y)- After this operation, we obtain 

A' 2 (x,yr E = + :h^(x)h^(y): a' 2 (x - y)%% , 
R'2(x,yr E = + :h^(x)h^(y): r' 2 (x - y)g^ ■ 

The C-number tensorial distributions a' 2 (x — y) 9 ^^ and r' 2 (x — y) 9 J^ u are given 
by linear combinations of products between positive/negative frequency parts 
of the Jordan-Pauli distributions (2.6) and carry four derivatives (each photon 



field in carries one derivative). Therefore, the basic distribution appearing 
here is 

" V) ■■= ^D^ix - y) ■ dffiD$\x - y) . (4.3) 
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Summing up all the contributions coming from —T£(x)Ti(y) after performing 
two photon contractions, we have 



+ 2 ^;l + + 



fip\pp 'PP ctp\pp 



(4.4) 



as a function of the D | -distributions. The distribution r^x — yfapiu/ nas ^ ne 
same form as above, but with D | instead of (because of the relation 

Dq (y — x) = —Dq(x — y)). The D , -distribution are easily evaluated in 
momentum space by means of Eq. (|2.6| ) : 



a.fi\\iv y r I 



■1 



(2vr) 4 

X 



ci 4 fc<5((p - A:) 2 ) 9(± (p° - A; )) ctyfc 2 ) e(±fc°) 
"b PaPf)kpk u — Pakf}k^k u — ppkak^kp + k a kf}k^k v , (4.5) 



(see App. 1 of |L3|]), because products of Jordan-Pauli distributions go over into 
convolutions of the corresponding Fourier transforms. Therefore, we see that 
we have to deal with integrals of the type 



I {±) (p)-/a/afS/al3p/a^ ■= / d*k 5((p - kf) 0(± (p° - k )) 5(k 2 ) Q(±k°) 

x [l, k a , k a kp, k a kpk^, kakpk^k^j . 



(4.6) 



These integrals have been calculated in detail in App. 2 of [ 13 1 and partially 
in [H|. We give here only the final results: 



7T 



i (±) (p) = ^e(p 2 )&(±P ) 



i {±) {pU 
i {±) (pU, 



PaP/3 



PaPpPfi 



vap)e(p 2 )®(±P ), 



(4.7) 
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Pa Vf3p + P/3 Vap + Pp Va/3 ) ®(P ) ®(±P 



I {±) {p)aPpu = ( PuPpPpPv - — ( + PaPfl T]pv + PaPp T]p u + p a p„ T]f3p, + 



+ PpPp Vau + PjiPv Vap + PpPv Vaj3 J + 

+ 48 ( + 7]a f MV ' 3u + Vauripu + VapViJ)) @(P 2 ) @(±P°) 
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From ([4.5D, (4.6) and (4.7), we obtain the relations for the l) 1 *,' -distributions 



£> (± > (p) 



1 



(2vr) 4 

— 7T 

48(2vr) 
-1 



■I- 

7 [2/^+/^] G(p 2 ) 8(±p°) 



(2vr)4 

— IT 



(4.8) 



7PJ 



8(2vr)' 



■p A e( P 2 )Q(±p° 



Therefore, from (gT2j) with G(p 2 )e(-p°) - G(f> 2 )9(+p°) = -6(p 2 ) sgn(p°) we 
find 



D 2 (x, V y SE =: h^(x)h^(y) : d 2 (x - yf SE 



a/3^tu ' 



(4.9) 



With Eqs. ( t4.4|) , (4.8) and after symmetrization in (a(3) «-> (/iz^) the C-number 
tensorial ^-distribution reads in momentum space 



K 2 7T 



960(2vr) 4 



- 16 PaPpPtiPv ~8p 2 {PaPpll^u + PnPuVa/3) + 

+ 12 p 2 (paPp.il/3u + PaPv-nPii, + PpPfiVav + PpPvVap) + 

- 12 p 4 (rj atl r]i3 U + rjavTipfj) 



+ 8p Vaprj^v 



&(p 2 )sgn(p°) 



(4.10) 



For simplicity, we use the shorthand notation 

^(P)S* := [-16, -8, +12, -12, +8] , 



(4.11) 



where T := k 2 7t/(960(27t) 4 ), and d(p) := @(p 2 )sgn(p°) (this scalar distribution 
is typical for a loop with massless particles) and the numerical coefficients always 
refer to the polynomial structure as in Eq. ( 4.10 ). The polynomial P{p) 9 ap E ^ u nas 
degree four in p, because of the four derivatives present on the two contracted 
photon lines. 



4.2 Singular Order, Distribution Splitting and Graviton Self- 
Energy Tensor 

According to the inductive construction of T 2 (x,y), Sec. || the next step is 
the splitting of D 2 (x,y) into a retarded part R 2 (x,y) and an advanced part 
A 2 (x, y). In this procedure the singular order of the distribution D 2 (x, y) plays 
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an essential role. From ( |4.1C| ) it follows that u^cfc) = 4, because of the presence 
of the polynomial of degree four in p. 

A general formula for the singular order of any distribution in linearized 
gravity coupled to photon fields can be found by considering in the n-th order 
of perturbation theory an arbitrary n-point distribution 



n h n A 



T n G (xi,... ,x n ) =:JJh(x fcj ) [jA(x mi ): t°(x lt ... ,x n ). (4.12) 

3=1 i=l 

This T° corresponds to a graph G with nh external graviton lines and n A ex- 
ternal photon lines. Then we state that the singular order of G is given by 

uj(G) <A-n h -n A - d + n. (4.13) 



Here, d is the number of derivatives on the external field operators in ( 4. 12] ). 



The explicit presence of the order of perturbation theory renders the theory 
'non-normalizable': the number of free undetermined, but finite normalization 
terms in Eq. ( |3,6[ ) increases with the order of perturbation theory, that is the 
theory has a weaker predictive power but it is still well-defined in the sense of 
UV finiteness. The proof of (|4.13Q has the same structure as in QED, Yang-Mills 
theories and pure quantum gravity, see therefore Q , || and ]l3[, r espectively. 

In the case of the graviton self-energy contribution, Eq. (|4.9| ), we obtain 
from fl4.13| ) the correct result ufy) = 4, being = 2, n A = 0, d = and 
n = 2. The singular order of a distribution remains unchanged after distribution 
splitting. 



Because of the decomposition (|4.1l[) , it suffices to split d(p), which has 



u){d) = 0, in order to obtain a retarded part f(p). The full retarded distribution 
is then given by r2(p) g a S g E = r(p). The ambiguity in the normaliza- 



tion appearing in the splitting formula (p.4|) will be discussed in Sec. 4.4. 



The splitting of the scalar distribution d(jp) was already carried out in ||], 



using a modification of the formula ( |3.5[) for the retarded part, see also [13|. We 
quote only the result for the analytic continuation to p 6 I 4 of the retarded 
part: 

f(rt °" = i log ( Z T 2 )' (414) 

where M > is a scale invariance breaking mass. Then 

R 2 (x,yy SE =:h^(x)h^{y): r 2 (x - y)%% , 

The distribution T2(x,y) is obtained from R,2(x,y) by subtracting R' 2 (x,y), 
Eq. (|4.2|). This subtraction affects only the scalar distribution. Since f'(p) = 



13 



-Q(p 2 ) 6(— p°), we obtain 



Hp) = f(p) an - r'ip) = ± log ( P ]'° ) . (4.16) 



2tt V M2 

Therefore, the 2-point distribution for the photon loop graviton self-energy reads 

\9SE 



T 2 (x, y y SE =:h^{x)h^{y): iU(x - y)f E 



p2- i0 \ (4.17) 



i n(p)^ = «3 [ — 16, -8, +12, -12, +8] log ( 

with 3 := 2^~T. The most important aspect of this result lies in the fact that 
the obtained T 2 (x, y) 5SB -distribution is UV finite without the introduction of 
counterterms ||, || and and cutoff- free. Note that M is a normalization 
constant and not a cutoff. 

4.3 Slavnov— Ward Identities from Perturbative Gauge Invari- 
ance 

We investigate the gauge properties of the graviton self-energy tensor. For 
simplicity let us denote by [A, B, C, E, F] the five numerical coefficients of Tl(x — 
y) a /3fiu given as in Eq. Q4.10D - 

First of all, T2(x,y) 9SE satisfies the condition of perturbative gauge in- 
variance to second order (|2.21 ). Using the infinitesimal operator gauge vari- 



ation ( |2.12| ), we obtain 

d Q T 2 (x,yY SE = + ^( :u"{x)hr{y): [ + b a ^ U(x - y) a ^ u ]) + 
+ %( :h^(x) U P(y): [ + V"" U(x - y) a ^] 



(4.18) 



because, after working out explicitly in momentum space, the necessary condi- 
tion for (pLilj ) to hold, 



Pa^(p)af3pu=0, (4.19) 



corresponds to the identity ( 2.27] ) and is satisfied by f[{p) a pn V of Eq. (14171 ) . In 



terms of the A, . . . ,F coefficients, (|4.19j ) implies 

C + E = 0, A-2B = 0, B -2E -2F = 0. (4.20) 

and these relations are satisfied by the polynomial in ( |4.17| ). 

We show now that these relations implies the gravitational Slavnov- Ward 
identities (SWI) [||, |nj. We construct the 2-point connected Green function 
with one photon loop 

G(P) 1 X» ■= W D F (p) fl{pT Spa b pafMU D E (p) , (4.21) 
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with the Feynman propagator Dq{p) 
SWI reads 

P a G(p) 



(2ir) (—p—iO) . The gravitational 



[2] 

a/3fj,u 



o, 



(4.22) 



namely the 2-point connected Green function is transversal. Eq. ( 4.22| ) implies 
the relations 



A 

C + E = 0, A-2B = , - + C - F = 0. 



(4.23) 



Analysis of (4.20) and (4.23) shows that the condition of perturbative gauge 
invariance is equivalent to the gravitational SWI. In addition, the photon loop is 
transversal also without the 6-tensors that come from the graviton propagators: 



because the equivalent relations 

C + E = 0, B + F = 0, A + B + 2C = 0. 



(4.24) 



(4.25) 



are fulfilled by ( 4.17| ). This property follows directly from our perturbative 
invariance condition (|2.27l ), because the trace of tl{p) a p^ v vanishes (see below). 

The last two properties of tl(p) a /3fMu concern its trace and are related to 
conformal transformations. Conformal invariance is manifest in the vanishing of 
the trace of the Maxwell energy-momentum tensor in four dimension. Basically, 
the photon loop graviton self-energy consists of a 'time-ordered' product of two 
such traceless photon energy- momentum tensors ( j2.17|) . Therefore, in addition 
to the identities ( |4.ig| ), it is expected that ti{p) a p^ v and G(p)^ are traceless. 
The condition 



7^11(^)^ = 0, 



(4.26) 



namely that the photon loop graviton self-energy tensor is traceless, implies the 
relations 



B + 2E + AF = , A + AB + AC = . 



(4.27) 



These are satisfied by ( [1.17 ). Note that ( 4.27] ) does not follow from the pertur- 
bative gauge invariance conditions (4.20). In addition, in our approach, also the 
trace of the 2-point connected Green function vanishes: 



because the equivalent relations 



— + 3B + 2C + 2E + AF 



, A + AB + AC = . 



(4.28) 



(4.29) 
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are satisfied by the coefficients of the photon loop self-energy tensor in Eq. ( [4.17] ) . 
This is in sharp contrast to the calculation in || , |ll]] , performed within the di- 
mensional regularization scheme. This scheme generates conformal trace anoma- 



lies [18], because the trace operation is not dimensional invariant. As a conse- 
quence, the extraction of the finite part from UV divergent quantities breaks 
this symmetry. 

Causal perturbation theory not only provides us with UV finite results but 
preserves the underlying symmetries of the theory such as the SWI and the 
vanishing of the trace of the self-energy tensor. Invariance under conformal 
transformations is broken by the presence of the mass scale M in the logarithm 



of Eq. (p?7D . If necessary, we could exploit the ambiguity in the normalization 



to restore it, see Sec. 4.4 



4.4 Freedom in the Normalization of the Graviton Self-Energy 

We still have to discuss the ambiguity in the splitting procedure, namely the 
appearance of undetermined local normalization terms. 

Having singular order four, the 2-point photon loop graviton self-energy 
contribution T2(x,y) 9SE admits a general normalization of the form 

N 2 (x, yf SE =: h^{x)h^{y) : i N(d x ,d y ) af}fiv 5 w (x - y) , 

n(pU, u = iv»lV + + #(p)&* > ( 4 - 3 °) 

where the odd terms are excluded by parity. N{p)^p is a polynomial in p of 
degree i with i = 0, 2, 4. We assume in addition that only scalar constants should 
be considered, because vector-valued or tensor-valued constants may endanger 
Lorentz covariance. Taking the symmetries of — y) a piii> into account, we 
make the following ansatz 

N (P)altiv = 3 [ + c l {VafiVpv + riavrip^) +c 2 rj 

^(P)%^ = H [ + C 3 (paPpV^ + P^Pur] a p) + 

+ C A [PaPuVflv + PaPuTIp^ + P(lP^ au + PpPuVafl) + ( 4 - 31 ) 
+ C 5 P 2 [VanVPv + VavVpii) + c 6 VapVfiv] , 

N(p)afa v = 3 [or, cb, as, cio, en] . 

ci, . . . , en are undetermined real numbers. The normalization polynomials have 
to fulfil the same symmetries as n a/ 3^, namely: (|4.2C| ), ( [4.25]) , fl4.27|) and (|4.29| ). 
These reduce the choice to 

*(P)SU = 0, N(p) ( Xv = 0, 
N(P) { X V = H C -f [ - 16, -8, +12, -12, +8] . (4.32) 

Since N(p) *g has the same structure as n Q /3/^> we absorb the undetermined 
parameter c\\ in the mass scale M appearing in Eq. (|4.17| ) through the rescaling 
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en = 41og(M 2 /Mg). The whole freedom in the normalization is reduced to the 
mass parameter Mq in the logarithm owing to the symmetries that the self- 
energy tensor and its normalization have to fulfil. 

This normalization automatically preserves graviton mass- and coupling 
constant-normalizations. If one sums up the infinite series with an increasing 
number of photon loop graviton self-energy insertions, the vanishing of N{p)^ 
implies that the graviton mass is not shifted by quantum corrections and the 
vanishing of N(p)*p implies that the coupling constant k is not changed by 
quantum corrections. 



As pointed out at the end of Sec. 4.3, we could exploit this normalization 



to restore conformal invariance. Through a rescaling of Mq we can compensate 
the variation of p 2 under such a transformation. The drawback lies in the 
arbitrariness of this operation. 

4.5 Corrections to the Newtonian Potential 



As pointed out in IS], [20], [21 and p^j , massless particle loop corrections to the 
graviton propagator leads to quantum corrections of the Newtonian potential 
between massive spinless bodies in the static non-relativistic limit. These can be 
appropriately defined by considering the whole set of diagrams in the scattering 
(piy>2 — > where ifi represents a scalar field of mass m», of the order k 4 (i.e. 

~ G 2 ). Then one isolates the non-local contributions. These lead to r~ 2 and 
r~ 3 corrections to the Newtonian potential V(r) = —Gmirri2r~ 1 in the static 
non-relativistic limit. For this purpose, the logarithm-dependent result ( 4.17| ) 
generates the correction 

^ =-G— ( 1 + — T7— ) , 4.33 
r \ c^ir 15 r z J 

as calculated in [O], where h and c are put back in the expression and the 
relevant length scale appears to be the Planck length l p = \J Gh/c 3 . 

Note that the mass scale M in Eq. (4.17) is irrelevant, because it contributes 



only to local terms ~ <5 (3) (a:). Therefore, the still remaining freedom in the nor- 
malization, namely the choice of M, is irrelevant to physical quantum corrections 
of the Newtonian potential. 

The correction in Eq. ( [4.33|) is only a partial one, because we have taken 
into account only the photon loop graviton self-energy contribution and not the 
complete set of diagrams of order k 4 contributing to these corrections, as, for 
example, the vertex correction or the double scattering. Therefore we cannot 
make any statement on the absolute sign and magnitude of the numerical factor 
in Eq. ( [4.331 ). To our knowledge, photon loop corrections to the exchanged 
graviton have not been considered yet. 
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5 Photon Self-Energy 

We now turn to the calculation of the photon self-energy. Since the main features 
of the causal scheme have already been presented, we can go on more speedily. 
Moreover this graph is not interesting from the point of view of gauge invariance: 
it is trivially gauge invariant due to d,Q A F a/3 = 0. 

In order to obtain the relevant A' 2 (x, y) and R' 2 (x, y), we carry out one gravi- 
ton and one field-strength contraction between T±(x) and T^(y) and obtain 

A' 2 (x,y) =:F^(x)F^(y): a' 2 (x - y)^ , (5.1) 

with 



k 2 

a' 2 {x - y) ay \ w ■= -r ( " V-yp - 7? w - r, ap + 



4 V 'hP^-\/ia 'hP^-\ap 'lap^-^n 

- 5 Vafl D ( J + 2r? 7a L>« + 2 Vw D^J (x - y) , (5.2) 

where we have as in (fl.5|) 



A c Ji(p) = ^yi l d4k6 ((p - k ?) ®(p° - k °) ■ (5-3) 



The expression in (|5.2| ) has to be antisymmetrized in a <-> 7 and in /i 

■I 



p, because these manifest antisymmetries in (|5.1|) get lost using the Z? (±) 



distributions. Therefore, with ( |4.6| ) we find in momentum space 

%l(p)ary\np = 48 ( 27r )4 [ + 2 PaP P rj 7l x ~ Zp-yPpVap ~ ^PaP^P + 2 PjPp,Vap + 

~ p\»Vap + P 2 Va f iV-yp] ©(p 2 ) ®(P°) • 

(5.4) 

The result for R' 2 (x, y) is analogous to that of A' 2 (x, y): the C-number tensor dis- 
tribution r' 2 (x — y) a -y\pp depends on the D , -distributions so that in momentum 
space the factor 0(— p°) is present instead of Q(p°). The causal /^-distribution 
is then given by 

D 2 (x,yr E =:F^{x)F^(y): d 2 {x - yf^ , 

foKw = t the same as in 0] (^ 2 ) S S^(P°) ■ (5-5) 

From Eq. ( [4.13 ) or from direct inspection, we ascertain that the singular or- 
der of the distribution is two. The distribution splitting is the same as in 
Sec. [l^, Eq. ( 4.16j ), because the loop consists of massless particles. Therefore, 
the T 2 (x, ^-distribution corresponding to the photon self-energy reads 

T 2 (x,yT SE =:F^(x)F^(y): ( - iU(x - yf^J , 
^(P)aZ P = 48^)5 t the Same aS in ©] l0 § ( P M ^ ) ■ (5 ' 6) 
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C(l) BM -perturbative gauge invariance to second order holds trivially: 

d QA T 2 (x,yr SE = 0. (5.7) 

The photon self-energy 2-point distribution can be written also in terms of the 
photon field A 11 : 

T 2 {x,yy SE =:A^x)^{y)y. ( - 4iH(s - yfJ^J 
=:A\x)A»{y): ( - i U(x - y) Jp ) , 

by disregarding non-contributing divergences (in the adiabatic limit g — ► 1 of 
Eq. (2.1)). Then T2(x,y) pSE has the same structure as the photon self-energy 
in QED Q and the 'reduced' photon self-energy tensor appearing Eq. ( |5.8[ ) is 

n(p) w = 4 P yn( P ) ai | W 

= 12W [p Wp " p nip] g V^^~ J 

and satisfies the Ward identity 

p'tl( P ) lp = (5.10) 

as in QED. 

As a last point, we discuss the normalization of T2(x,y) pSE in the second 



line of Eq. (|5.8| ) with singular order four from (5.£). Actually, we should con- 
sider the normalization of Eq. (5.6) with singular order two, but owing to the 
relation ( |5.9[ ), the two analyses lead to the same conclusion. The normalization 
terms have the form 

N 2 (x, yf SE =:A^{x)A"{y) N{d x ,d y ) lp 5 w (x - y)) , 

N(p) lp = iV»<°> + N(p)V + N(p)% . (5.11) 

Due to the same reasons as pointed out in Sec. [O], the local normalization can 
be expressed in momentum space through the polynomials 

N(p)% = * [ Cl r, 7P ] , 

Nip) 1 * = [c2 PlPp + c 3 ^ pP 2 ] , (5.12) 

Here, c±, . . . ,05 are undetermined real numbers and ^ := K 2 7r/12(27r) 5 . Requir- 
ing (|5.10| ) for N(p)typ, i = 0,2,4, we reduce the freedom in the choice of the 
normalization parameters in the polynomials to 

JV(p)$ = 0, N(p)% = C2*[p lPp -riwP 2 }, (5-13) 
#(p)$ = c 4 [p 2 PyPp - W 4 ] . (5.14) 
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The normalized photon self-energy tensor then reads 



n(p) 



N 

IP 



Hp 2 )n := * 



p 4 log 



fh/p) U (P 2 



with 



-p 2 - iO 
M 2 



+ c 2 p 2 + c 4 p 4 



(5.15) 



In order to fix the remaining free parameters c 2 and C4, we consider the total 
photon propagator in momentum space, defined as the sum of the free photon 
Feynman propagator and an increasing number of self-energy insertions [fjj 



+V,u D F (p) + D F (p) n(p)* D F (p) + 

+ D F {p) TL(p)» x D F (p) fl(p)t D F (p) + ... 

b F ( P ) [ Vllu + n( P )» x D(p)^], 



(5.16) 



with fi(p)£, := (2vr) 4 fl(p)^. Multiplying (fTl6|) with Dg(p)- 1 and with 77^ we 
obtain 

[V ^(p)" 1 " n(p)^ A ] D(p)g* = ry w . 
The inverse of the total photon propagator is then 

(Dip)^ = (2vr) 2 [rj^-p 2 - <0) - (2vr) 2 n(p)*] . 
Therefore, inverting the expression above, the total propagator reads 



D{p) 



tot 



(2vr)< 



-p 2 - io - no? 2 ), 



, P^P" T?l 2 
+ —2-^0 



(5.17) 
(5.18) 

(5.19) 



The explicit form of the function F(p 2 ) is not important, because the last term 
vanishes between transversal photon operators. Photon mass- and coupling 
constant normalization 



p 2 =0 



and 



p 2 =0 



0. 



(5.20) 



yield c 2 = 0. The last parameter C4 is non-essential, because it only shifts the 
mass scale M. Therefore, all local finite ambiguities in ( |5.11| ) can be reduced to 
a single unknown normalization parameter. 

This concludes our discussion of the loop graphs to second order in causal 
perturbation theory. 
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